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Transverse momentum spectra of pions, protons and antiprotons in Au+Au collisions at interme- 
diate RHIC energy of -y/ijviv = 62 GeV are studied in a model that includes both quark coalescence 
from the dense partonic matter and fragmentation of the quenched perturbative minijet partons. 
The resulting baryon to meson ratio at intermediate transverse momenta is predicted to be larger 
than that seen in experiments at higher center of mass energies. 
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I. INTRODUCTION 



Recently, there is a renewed interest in using the quark 
coalescence or recombination model to study hadron pro- 
duction in ultra-relativistic heavy ion reactions Q, |^ 
^ IE IE IE Hi 1^ ' These studies were largely triggered by 
two surprising observations in measured hadron spectra 
from Au-I-Au collision at ^/sNN = 200 GeV at the Rel- 
ativistic Heavy Ion Collider (RHIC) QUE EE El- The 
first was the nearly similar yield of protons and pions 
at intermediate transverse momentum (2 GeV < pt < 
5 GeV), which is at variance with the expectation from 
the fragmentation of minijets produced from initial hard 
processes. The second was the so-called "quark number 
scaling" of hadron elliptic flow V2h{pT)i i-e., the trans- 
verse momentum dependence becomes universal if both 
V2h and transverse momentum px are divided by the 
number of constituent quarks in the hadron. Both phe- 
nomena have a simple explanation if hadronization of the 
thermally equilibrated and collectively flowing partonic 
matter formed in these collisions proceeds through the 
coalescence of quarks of constituent masses |E S IE 13] ■ 

The description of hadronization of the quark-gluon 
plasma formed in heavy-ion collision at relativistic ener- 
gies by quark coalescence was first introduced in models 
like ALCOR 13] and MICOR Jl]. Emphases of these 
earlier studies were on particle yields and their ratios. 
Recent studies have been more concerned with observ- 
ables that are related to collective dynamics and produc- 
tion of hadrons with relatively large transverse momenta. 
Furthermore, effects of minijet partons from initial hard 
processes have also been included through their indepen- 
dent fragmentation as well as coalescence with soft par- 
tons in the quark-gluon plasma jE IE] ■ The interplay be- 
tween quark coalescence and minijet fragmentation has 
been found to be essential for understanding measured in- 
clusive hadron spectra at moderate and high transverse 
momenta. If confirmed by more exclusive data from ex- 
periments and by further theoretical studies, hadroniza- 



tion via quark coalescence then represents a new probe of 
the quark-gluon plasma formed in relativistic heavy ion 
collisions. It is thus of great interest to study the onset 
of the coalescence mechanism as the collision energy is 
varied. Also, comparison of results from the coalescence 
model for heavy ion collisions at different energies with 
upcoming experimental data will allow for a better un- 
derstanding of its predictive power. In the present paper, 
such a study is carried out by using the same coalescence 
model as applied previously to higher energy collisions 
to make predictions for the proton, antiproton, and pion 
spectra at a lower collision energy of ^/sWn = 62 GeV. 

This paper is organized as follows. In Sect J^we briefly 
review the coalescence model used in Ref. 0. How the 
parameters in the model are determined for Au-|-Au col- 
lisions at y/SNN = 62 GeV is described in Sect. IIIll We 
then discuss in Sect. llVl the fragmentation contribution 
from a pQCD approach that takes into account both 
the Cronin effect and the radiative in-medium energy 
loss EE]- In Sect. O results on pion, proton and an- 
tiproton spectra and their ratios are presented. Finally, 
a summary is given in Sect. IVII 



II. THE COALESCENCE MODEL 

The coalescence model describes the dynamical process 
of converting quarks and antiquarks into hadron bound 
states in the presence of partonic matter. In general, the 
process of coalescence would involve the emission of a 
third particle or, more generally, off-shell effects to guar- 
antee energy conservation and color neutrality. However, 
under the assumption of a fast process in which the bind- 
ing interaction is turned on suddenly, the probability of 
forming the bound state is simply given by the overlap of 
the quark and antiquark distribution functions with the 
Wigner function of the formed hadron. In this case, the 
transverse momentum spectrum of a hadron consisting 
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of n (anti-)quarks can be written as 

In the above, da denotes an element of a space-like hyper- 
surface, and gn is the statistical factor for forming a col- 
orless hadron from colored quarks and antiquarks of spin 
1/2. For hadrons considered here, i.e., it, p, K* , p, and A, 
the statistical factors are, respectively, = g^j = 1/36, 
9p = gK' = 1/12, gp = 1/108, and = 1/54. The 
n-quark phase space distribution function has been ap- 
proximated in Eq. ^ by the product of the single quark 
distribution function fq{x,p) that is normalized to the 
quark number in the partonic matter. 

The Wigner function fH{xi;pi) in Eq. describes 
the spatial and momentum distributions of quarks in a 
hadron. For mesons, we have taken the Wigner functions 
to be a sphere of radius A^, in coordinate space and Ap — 
in momentum space, i.e.. 
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(2) 



This represents the simplest ansatz for the Wigner func- 
tion, i.e, it gives a probability of one if the relative covari- 
ant distances between two particles in space and momen- 
tum are smaller than the radius parameters. The Wigner 
functions for baryons can be similarly expressed in terms 
of appropriate relative coordinates (see Ref. Q for de- 
tail). Although hadron wave functions are generally not 
known in detail, the particular shape of the wave function 
used in our model does not affect significantly the result 
on hadron transverse momentum spectra, unless very ex- 
treme choices are taken. As in our previous work |2j, we 
take Ap = 0.24 GeV for mesons and Ap = 0.36 GeV for 
baryons in the present study. 



III. FIREBALL PARAMETERS 

One of the main ingredients in our model is the par- 
ton distribution functions at hadronization. They can 
in principle be determined from a d yna mic approach 
such as the transport model E mill El llfl. To 
focus on the hadronization process itself, we assume in- 
stead that at hadronization the partonic matter consists 
of both a thermally and chemically equilibrated quark- 
gluon plasma and quenched minijet partons from initial 
hard processes. Specifically, the single quark distribu- 
tion function is taken to consist of both a thermal and 
a minijet component, separated by the transverse mo- 
mentum po ~ 2 GeV. In this section we discuss how 
the thermal component is determined from the mea- 
sured total transverse energy as in our previous work 



at ^/snn = 200 GeV [^, while the minijet one will be 
discussed in the next section together with the minijet 
fragmentation process. 

Under the assumption that the quark-giuon plasma is 
thermally equilibrated and collectively flowing, the deter- 
mination of its parton distributions requires the knowl- 
edge of the radial flow and its total transverse energy. 
A blast wave analysis of the experimental data shows 
that the collective radial flow differs very little between 
the highest SPS and RHIC energies. We can there- 
fore assume the same radial flow for Au-|-Au collisions 
at Y^sjvw = 62 GeV as at higher RHIC energies, i.e., 
Pmax — 0.5 or a slightly smaller one of about 0.45 and 
the usual linear radial profile /3 = Pmax f / R, where R 
is the transverse size of the fireball. We have checked 
that this uncertainty does not have a significant effect on 
the transverse momentum dependence of the p/tt ratio, 
which is the main observable we want to study in this 
paper. 

As to the total transverse energy, we are not aware at 
the moment of any measurement for Au-|-Au collisions 
at y^SNN = 62 GeV. It has been, however, noticed that 
the ratio of the transverse energy to the total number 
of charged particles in relativistic heavy ion collisions is 
nearly constant in the energy range 20 GeV < ^Sjvjv < 
200 GeV, and available experimental data [IJ can be 
parameterized by 
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0.88 -f 1.44- 10-3(7i^- 130) GeV. (3) 



Also, the charged particle multiplicity normalized to the 
number of participants has a smooth logarithmic behav- 
ior that can be parameterized as 
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Using the Glauber model to estimate the number of par- 
ticipants, we have Npart = 330 for the — 10% central- 
ity bin. From Eq. Q), the charged particle multiplicity 
is then dNch/drj — 480, corresponding to a multiplicity 
per unit of rapidity of dNch/dy = 1.2 dNch/drj — 580, 
with drj/dy — 1.2 from the relation between rapidity and 
pseudorapidity. Eq. ^ thus leads to a transverse en- 
ergy dEr /dy ~ 460 GeV for central Au-I- Au collisions at 
= 62 GeV. 

What remains to be determined is the relative abun- 
dance of quarks and antiquarks, which is essentially given 
by the baryon chemical potential. This can in prin- 
ciple also be extracted from experimental data on the 
p/p ratio via fig = — T/6 ln(A'p/A'p). For collisions at 
y/sNN — 200 GeV, we already had an experimental value 
of p/p ~ 0.7, which corresponds to a quark chemical po- 
tential /i, = 11 MeV |3|. There is no measurement of 
the p/p ratio yet at ^/sNN = 62 GeV, but it can be esti- 
mated from a parameterization of the chemical potential 
in the energy range of 10 GeV < y/sJfN < 200 GeV, 
based on available data from AGS to RHIC. The result- 
ing quark chemical potential is fig — 23 MeV, leading 
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to an antiproton to proton ratio p/p= 0.44. Because of 
strangeness conservation, numbers of strange and anti- 
strange quarks are equal, i.e., zero strangeness chemical 
potential. The relative abundance of strange to light 
quark is then simply determined by their constituent 
quark masses that are taken to be mu.d — 300 MeV 
and TTis = 475 MeV. Assuming a critical temperature 
of Tc = 170 MeV, the partonic matter at hadroniza- 
tion is found to have Nu+d = 285, -/V^+j = 218, and 
Ng — Ng — 79 quarks and antiquarks distributed in a 
volume V=555 fm"^ in order to have the estimated trans- 
verse energy of 460 GeV. This leads to a critical energy 
density of about 0.8 GeV/fm'^, similar to that from the 
lattice QCD calculations 



IV. QUENCHED PERTURBATIVE QCD 
PARTONS AND HADRONS 

Partons with transverse momenta pr > Po (few GeV) 
are mainly produced from initial hard collisions between 
nucleons. Their distributions can be determined from 
the pQCD calculation that includes the Cronin effect and 
radiative energy loss as described below. 

In "elementary" hadron-hadron {N + N) collisions, 
high pt minijet parton production can be evaluated us- 
ing the perturbative QCD factorization approach [23l |. 
This leads to a differential cross section that is expressed 
as a convolution of the measured parton distribution 
functions 24] fa/N{xa,Q\) for the interacting partons 
(q — a, b) with the elementary parton-parton cross sec- 
tions da^"^^"'^'^ / di 
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X fO^Ta)fO^Tb)fa/p{Xa, Ql)fb/p{xb, Ol) 
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TT dt 



-5{s + u + t). 



(5) 



In the above, Xa and Xb are the initial momentum frac- 
tions carried by the interacting partons. For nucleon- 
nucleon collisions, we use phenomenological smearing 
(k^) — 1.8 GeV^/c^ via a normalized Gaussian distribu- 
tion function /(kx) to mimic next-to- leading order and 
soft gluon resummation effects. The determination of 
Kn lo in the partonic cross sections in the absence of ini- 
tial and final state multiple scattering is done as in [26l |. 

In the presence of cold and hot nuclear matter the mul- 
tiple elastic, inelastic and coherent scattering of the in- 
coming and outgoing partons modify the perturbative 
cross sections Eq. ((SJ. The calculated corrections can 
be systematically incorporated in the pQCD factoriza- 
tion approach 27] . For example, transverse momentum 
diffusion in p + A collisions leads to the Cronin effect 
|26j . In A -\- A reactions the dominant nuclear effect is 
the energy loss of energetic jets propagating through the 
quark-gluon plasma and the dominant energy loss mech- 
anism is via medium-induced gluon bremsstrahlung. To 



first order in the mean number of scatterings, the radia- 
tive spectrum and the average energ y lo ss, calculated in 
the Gyulassy-Levai-Vitev approach [23 for the case of 
Bjorken expansion reads 



dN. 



9 

ind. 



dui 



Ci 



97ra2 1 dNP° 



4 A 



dy 



L 



UJ 



— In ^ + • • • , « t?L 2 



(6) 



1 



Qual 1 dN^" 



In- 



2E 



A^ 
3 



dy 



L 



(7) 



where Cr = 4/3 (3) for quark (gluon) jets, A±^ is the 
transverse size of the interaction region, and L is the 
distance traversed by the jets. We note that for the pur- 
pose of calculating the attenuation of the single inclusive 
hadronic spectra 1+lD expansion approximates well the 
3-l-lD result. 

Eqs. © and illustrate qualitatively the behavior of 
the energy loss as a function of the system size, density 
and partonic energy. Quantitative jet quenching analysis 
requires numerical evaluation of the the radiative spectra 
and the mean energy loss as in Refs. In the Pois- 

son approximation of independent gluon emission, the 
probability P(e) for fractional energy loss e = X^i^j/^ 
due to multiple bremsstrahlung can be generated from 
the radiative spectrum as in Ref.js^]. Along the lines 
of Rcf. we choose a suppression factor that corre- 
sponds to an initial parton density dNparton/ dy = 725. 
The resulting quenched minijet distributions are used for 
the parton distributions at pT > Po in the coalescence 
process. 

Since the probability for a minijet parton to coalesce 
with other partons to form a hadron is quite small, in- 
dependent fragmentation remains the dominant mecha- 
nism for hadronization of minijet s 0|. This contribution 
is obtained by introducing the fragmentation function 
Dh/c{z, Qc) for minijet parton c into a hadron of flavor 
/i, i.e., 



dNhad 

dyd"^ Pt 
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dNjet{PT/zc) Dh/c{zc,Ql) 



where Zc — Ph/Pc is the momentum fraction carried by 
the observed hadron and Qc = Ph/zc is the momentum 
scale for the fragmentation process. It is this perturba- 
tively computed hadron spectrum that we refer as frag- 
mentation component of the hadron spectrum. 

For pions, we use the fragmentation functions of 
Ref. m that are determined from e+e , ep, and pp col- 
lisions. For pp collisions at ^/SNN — 62 GeV, where both 
the Cronin and radiative energy loss effects are absent, 
the predicted pion spectrum reproduces very well the ex- 
perimental data at > 3 - 4 GeV but overestimates that 
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at = 2 GeV by about a factor of 2. To also reproduce 
the pion yield at pt — 2 GeV, we have thus modified the 
minijet parton distribution at low px- For baryons such 
as p and p, we use a string-inspired parameterization of 
their fragmentation functions as in Refs. |i32i. .33]. The 
resulting p and p spectra are, nevertheless, sensitive to 
the quark and antiquark distributions in the nucleon. 

We note that application of the energy loss calcula- 
tions alone for neutral pions at ^/s =62 GeV in Au+Au 
collisions at RHIC [31 yields a factor 2 — 3 suppression 
with only a moderate dependence on px- The medium- 
induced non-Abelian bremsstrahlung thus plays an im- 
portant role in modifying the relative importance of the 
pQCD fragmentation mechanism when compared to the 
low pt non-pcrturbative approaches such as string and 
baryon junction mechanisms (33i. i34| . relativistic hydro- 
dynamics 1^ ) a-iid quark coalescence ^ IE ^| • 



V. PARTICLE SPECTRA AND RATIOS 

To obtain hadron transverse momentum spectra, we 
evaluate the multidimensional integral in Eq. ^ by the 
Monte Carlo method using test particles that are dis- 
tributed uniformly in the fireball according to the dis- 
tribution function and geometry described in Sects. 
and II VI Details of this method can be found in Ref. Q. 
As in our studies for Au-f Au collision at center of mass 
energy ^/snn — 200 GeV, we include quark coalescence 
into stable hadrons such as tt, p, and p as well as unsta- 
ble resonances such as p, uj, A, and K*. Decays of the 
resonances are included, and they modify the spectra of 
stable hadrons mainly at < 2 GeV, leading to a better 
agreement with experimental data at ^/snn = 200 GeV. 

In Fig.^ we show predicted pr spectra of pions (solid 
line), protons (dash-dotted line) and antiprotons (dashed 
line) including both coalescence and fragmentation con- 
tributions for Au-|-Au collisions at y^SNN = 62 GeV. For 
comparison, we also show the experimental pion (filled 
circles) and proton (open squares) data at = 
200 GeV from the PHENIX collaboration Is^. At low 
Pt, both pion and proton spectra are similar at the two 
energies. This is not surprising as the total charged par- 
ticle multiplicity, which is dominated by low pT parti- 
cles, is only reduced by about 20% as the collision en- 
ergy decreases from 200 to 62 GeV. At higher pT the 
pion spectrum depends mostly on the fragmentation of 
quenched minijets. With respect to the case at 200 GeV, 
the quenching mechanism is somewhat less effective due 
to a lower value of the initial gluon density (ISj . Prelimi- 
nary data from the PHOBOS experiment on the suppres- 
sion factor for charged particles has, indeed, confirmed 
such a trend 38] . On the other hand, less minijet par- 
tons are initially produced at ^/snn — 62 GeV than 
at y^s^vTV = 200 GeV and their transverse momentum 
distribution is also steeper, which amplifies the quench- 
ing effect at high px- The interplay between the non- 
perturbative and perturbative yields results in a tt spec- 
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FIG. 1: (Color online) Hadron transverse momentum spectra 
from central Au-I-Au collisions at ^/snn ~ 62 GeV. Predicted 
results are for pions (solid line), protons (dash-dotted line), 
and antiprotons (dashed line). Experimental 7r° (filled cir- 
cles) and proton (open squares) data at ^sjvjv = 200 GeV 
I39II are shown for reference. Inset: The ratio r{pT) of 
the contribution from quark coalescence to that from both 
quark coalescence and minjet fragmentation for tt (solid line), 
p (dash-dotted line), and p (dashed lixne). 



trum that is similar to the one at ^/snn — 200 GeV 
around px — 2.5 GeV but falls down more quickly as 
shown in Fig. ^ 

Fig. n further shows that although pions dominate at 
low Pt, more protons are produced at intermediate pt 
up to 5.5 GeV. The antiproton yield also becomes com- 
parable to the pion yield at intermediate transverse mo- 
menta. Since the minijet fragmentation gives the pro- 
ton and antiproton spectra that are about a factor of 4 
lower than the pion spectrum, the large proton and an- 
tiproton yields at intermediate transverse momenta are 
mainly due to contributions from the quark coalescence. 
This is further illustrated in the inset of Fig. ^ where the 
ratio r(j>T) of the contribution from quark coalescence to 
that from both quark coalescence and minijet fragmen- 
tation is shown for tt (solid line), p (dash-dotted line), 
and p (dashed line). It is seen that the quark coalescence 
contribution to proton production at intermediate trans- 
verse momenta overtakes the suppressed proton spectrum 
from minijet fragmentation up to px < 6.5 GeV, while 
for pions it is dominant only up to < 4 GeV. 

The predicted p/n ratio in central Au-|-Au collisions 
at ^SNN — 62 GeV is shown in Fig. [3 and it reaches 
a value of about 1.5 at ~ 4 GeV. For compari- 
son, the p/tt ratio of about 1 was measured in colli- 
sions at y^SNN = 200 GeV as shown by open squares 
for the experimental data [s^. The large value at lower 
collision energy is due to the larger difference between 
the slopes of the spectra from fragmentation and coa- 
lescence in collisions at ^/sNN = 62 GeV than that at 
y^SNN = 200 GeV. Also shown in Fig. |21is the p/n ratio 
at collision energy of y^SNN — 62 GeV. It can reach a 
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FIG. 2; (Color online) Same as Fig.Qfor the ratios of proton 
(solid line) and antiproton (dashed line) to pion spectra. 

maximum value of about 0.7, which is close to the value 
at y/sNN = 200 GeV, shown by filled squares for the ex- 
perimental data, but at slightly larger pT, as a result of 
the steeper tt spectrum from the fragmentation contri- 
bution. Because of the difference in the baryon chemical 
potentials, a larger difference between p/ir and p/tt ratios 
is seen at ^/snn — 62 GeV than at 200 GeV. 



at — 4 GeV, and a proton spectrum that is instead 
still dominated by contributions from the quark coales- 
cence up to Pt — 6.5 GeV. Experimental confirmation of 
this prediction will provide a stronger evidence for quark 
coalescence as a likely non-perturbative mechanism for 
hadronization, especially for protons (and in general the 
baryons) at transverse momenta up to about 6 GeV. This 
in turn will also facilitate the determination of the charac- 
teristic energy at which the quark-gluon plasma is formed 
in heavy-ion reactions. 

Our results depend on the properties of the bulk par- 
tonic matter at the moment of hadronization and the 
amount of minijet quenching. In the present study, 
they are determined from extrapolations based on avail- 
able data at ^snn — 17 GeV from SPS as well as 
= 130 and 200 GeV from RHIC. With upcoming 
experimental data on the total transverse energy, baryon 
chemical potential, and minijet suppression factor in col- 
lisions at y^s]vjv = 62 GeV, more reliable predictions on 
hadron spectra and their ratios can be obtained. Our 
work, however, represents an important first step in in- 
vestigating the interplay between the non-perturbative 
and perturbative hadron production mechanisms at in- 
termediate RHIC energies. 



VI. SUMMARY 

Using a hadronization model based on quark coales- 
cence and minijet fragmentation, we predicted the trans- 
verse momentum spectra of pions, protons and antipro- 
tons for central Au-|-Au collisions at ^/snn = 62 GeV. 
We found a larger enhancement of the p/tt ratio with 
respect to the one seen in nuclear collisions at ^^snn = 
200 GeV. Our result is mainly due to a steeper pion spec- 
trum, that is dominated by minijet fragmentation already 
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